Calculation of Cu/Ta interface electron transmission and effect on conductivity in 

nanoscale interconnect technology 
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Resistivity augmentation in nanoscale metal interconnects is a performance limiting factor in 
integrated circuits. Here we present calculations of electron scattering and transmission at the 
interface between Cu interconnects and their barrier layers, in this case Ta. We also present a 
semiclassical model to predict the technological impact of this scattering and find that a barrier 
layer can significantly decrease conductivity, consistent with previously published measurements. 



As the minimum feature size in transistors continues 
to shrink, nanoscale metal wires with thickness <32 nm 
will soon be needed to interconnect transistors in inte- 
grated circuits. However, measurements show nanoscale 
metallic wires have substantially higher resistivity than 
bulk metals [l|, Q , leading to performance-limiting inter- 
connect delays and power dissipation Q . Scattering from 
rough wire surfaces, interfaces with liner layers, and grain 
boundaries are believed to be the causes of this conduc- 
tivity degradation but microscopic understanding of 
these effects and quantitative predictions of their magni- 
tude have been limited. In recent works, we have consid- 
ered scattering from surfaces Q and grain boundaries [j| , 
and in this paper, we present calculations of scattering 
and transmission at the interface between Cu intercon- 
nects and their liner layers. 

Nanoscale Cu wires in integrated circuits are sur- 
rounded on three sides by a liner layer, typically made of 
Ta, TiN, or Ti. Also known as a barrier/adhesion/seed 
layer, this layer acts as a seed for the deposition of Cu, 
improves adhesion of the Cu to the sidewalk, and pre- 
vents diffusion of Cu into the dielectric. To inhibit the 
migration of Cu atoms, the liner layer is made of a re- 
fractory metal with small crystallites; this same prop- 
erty gives it poor conductivity despite being metallic. In 
this paper, we consider for the liner layer the most com- 
mon form of Ta when deposited epitaxially on Cu, (3— Ta 
0Ji which has a bulk resistivity of (200 ± 20) //f2- 
cm falsi, III, compared to 1.7 fj.fl- cm for Cu. 

has a tetragonal lattice with dimensions 
a = 10.2 A and c = 5.3 A. In the [002] direction (cor- 
responding to the c lattice parameter), it contains four 
equally spaced layers, which consist of either four regu- 
larly spaced Ta atoms or 11 atoms arranged in a pseudo- 
hexagonal pattern. The (3— Ta unit cell contains a total 
of 30 atoms. Simulations [l(| and measurements Q indi- 



cate that for Cu deposited on a (002) (3— Ta liner layer, 
the /3— Ta matches heteroepitaxially to Cu (111) surfaces 
with a relative strain in each a direction of approximately 
7%, causing the pseudo-hexagons in Ta to match those 
in the Cu (111) planes. 

We have set up such an epitaxial simulation cell for re- 
laxation with the VASP density functional theory (DFT) 
code [ill ] . Our interface contains a total of 56 Ta and 
64 Cu atoms with periodic boundary conditions in the 
plane of the interface (Figure QJ. These boundary con- 
ditions are made realistic by the epitaxial relationship, 
which should impose rough periodicity over actual in- 
terfaces. As shown in Fig. [TJ our Ta region terminates 
with an 11-atom Ta layer, interfacing with a partial Cu 
layer with four atoms, followed by four full (15-atom) Cu 
(111) layers. We considered a range of interface struc- 
tures, and this gave the lowest energy. We relaxed this 
system within DFT, keeping the bottom Ta layer fixed 
and allowing the top Cu layer to adjust only in the z di- 
mension. This allowed the interface z— spacing to relax. 
The variation in z— position among Cu atoms within the 
topmost layer after relaxation was very small compared 
to the layer spacing. 

To simulate transmission across a Cu/Ta interface, we 
replaced z— coordinates of the topmost Cu layer by their 
average, leaving ideal Cu and ideal strained (3— Ta layers 
bounding the interface region. We then set up a trans- 
mission simulation across the interface, using the relaxed 
interface of Fig. Q] as a scattering region and matching the 
perfect Cu and Ta layers to electrodes made of ideal lat- 
tices of the respective materials. We used the code Atom- 
istix 1121 ] to perform Non-Equilibrium Green's Function 



Method [ill simulations of dynamical transmission across 
interfaces, as in our previous work 

In this formulation, transmission probability 



T = T/M, 



(1) 
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where T is the total transmission (conductance in units 
of 2e 2 /h), and M is the Sharvin conductance in units 
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FIG. 1: Relaxed Cu/Ta interface, as used in transport sim- 
ulation. Top eight monolayers (darker, blue online) are Cu, 
bottom (lighter, brown online) layers are Ta. This structure 
simulates a Cu wire with Ta liner layer. The wire axis is 
parallel to the interface (left to right, as shown). The trans- 
port simulations use infinite Ta and Cu leads at the top and 
bottom of the system shown, and calculate the probability of 
electrons to transmit across the interface from Cu to Ta and 
Ta to Cu. 



of 2e 2 /h, which is equal to the number of current- 
carrying modes (l3j . Note that conduction electron den- 
sity, and therefore M, is different in the two materi- 
als. We obtained M for (002) strained /3-Ta and (111) 
Cu by separate transmission simulations, finding a ra- 
tio Mcu/Mth — 1-8. However, reciprocity requires that 
T is the same for transmission in either direction across 
the interface so the net current is zero across the 
unbiased interface after the vacuum levels have adjusted 
to equilibrate the Fermi levels We calculated T for 
transmission from Ta to Cu only, using the M's to find 
an overall transmission probability T = 0.22 for electrons 
originating in Cu and 0.39 for those originating in Ta. 

Let us consider the effect on conduction of transmission 
across the interface. Although the liner layer represents 
a parallel conductance to the Cu wire, the high resis- 
tivity of p— Ta makes this additional amount negligible 
14j . Instead, we consider the more important effect of 
electrons entering Ta and rapidly losing their net momen- 
tum before returning to Cu. If electrons encounter the 
interface, they have probability T of transmitting into Ta 
where they will lose their momentum. In steady state, 
for every electron leaving Cu there is on average one en- 
tering from Ta. Given the short mean free path in [3— Ta, 
the expected value of drift velocity for returning electrons 



is near zero. Effectively, electrons encounter the surface 
and lose their net momentum. Thus instead of paral- 
lel conductances, a more appropriate model for the Ta 
barrier layer is added surface scattering in a Cu wire. 

We can estimate the effect on conduction in Cu using a 
simple model with the same form as the surface scatter- 
ing models of Fuch s fTEj and Sondheimer [l(| . Here we 
follow Sondheimer [16| to find the semiclassical surface 
encounter rate (but we interpret the momentum loss as 
being due to electrons returning to Cu from Ta, rather 
than scattering at the surface). We assume only a frac- 
tion / of total carriers are affected by the surfaces, and 
as long as thickness is not too small compared to a mean 
free path A, only carriers within ~ A of the surface should 
be affected. The fraction becomes: 



/ = 



CPX 



where A and P are the wire cross section and perimeter, 
and we introduce a dimensionless constant C. Compar- 
ison with geometry-specific calculations gives C = 3/16 
[l6j]. The associated resistivity augmentation over bulk 
is given by 



P i 3AP ^ 
Pb 16^4 



(2) 



for arbitrary cross section, where 1 — p is the proportion of 
electrons that lose their momentum on interaction with 
the surface. In our case, this is the proportion of electrons 
transmitting into the barrier layer and losing their net 
drift velocity before returning to Cu. 

The assumption of zero drift velocity for returning elec- 
trons is valid when the Ta thickness is large compared 
with the bulk mean free path Ax a in P~ Ta, which we 
estimate as follows. The product of Sharvin (ballistic) 
conductance and bulk resistivity gives a length scale of 
order the bulk mean free path [1 31 ] , a relation that holds 
to within 40% for Cu[lJ. For strained (002) /3-Ta, our 
calculation of Mr a combined with bulk resistivity mea- 
surements suggest Axa ~ 1 nm (compared to A = 39 nm 
inCu). 

We now combine Eq. @ with our calculated transmis- 
sion probability to estimate the effect on conductivity 
from scattering in the barrier layer. Treating transmis- 
sion into Ta as equivalent to a diffuse scattering event, 
the effective (1 — p) for a wire would be given by 

1-JW=(1-P.) ^^ + {1-Pi + T} J, (3) 

where P is the total wire perimeter as in Eq. @, Pi is 
the perimeter interfacing with the liner layer, 1 — pi is the 
diffuse reflection probability from the Cu/Ta interface, 
and 1 — p s = 0.04 is our previously calculated diffuse 
probability for a rough Cu surface |4(. Here we will use 
Pi « p s as a crude approximation. Eqs. (H])-© give a 
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resistivity augmentation of about 13% for a 45 nm Cu 
square wire surrounded by Ta of more than about 2 nm 
on three sides, and 4% for a 45 nm Cu film with Ta on 
one surface. 

We can also easily extend this analysis to consider thin- 
ner liner layers of thickness txa ■ The cumulative scatter- 
ing probability is Poisson: 

P(s) = l-e- s ' x , 

where s is distance traveled and A is mean free path. 
Then the probability that an electron leaving Ta scat- 
tered since coming from Cu is given in terms of the trans- 
mission probability from Ta, Ty a = T/Mxai by: 
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(Here n is the number of internal reflections within Ta, 
which each occur with probability 1 — Ira)- So we find 



1 - Pi 



Pliaer(*Ta) = (1 - Ps) 

1 p — 2tTa/ATa 
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which interpolates smoothly between our surface 
scattering-only result and Eq. ([3]). Here Tq u — T/Mc u 
is the transmission probability from Cu. 

The form of Eq. (j4]) agrees qualitatively with the re- 
sults of Rossnagel and Kuan 14|, who measured sheet 
resistance of a 45 nm Cu film with 0-5 nm of Ta de- 
posited on it (note they use a film and deposit Ta on 
only one surface). They find resistivity augmentation 
rapidly increasing by 10% for ~2 nm of Ta, and leveling 
off for thicker layers. In Figure we show their data, to- 
gether with fits to Eq. using both our calculated value 
of T Cu = 0.22 and to the best-fit value of T Cu = 0.54. 
These two fits use Ax a ~ 1-2 nm, and Ax a ~ 1-8 nm re- 
spectively, for the mean free path in /3-Ta, in agreement 
with our expectation of ~1 nm. Although our calculation 
of Tcu fails to reproduce quantitatively the experimental 
curve in Fig. [2] the figure shows both that the form of ((4]) 
agrees with experimental data, and that we can predict 
behavior qualitatively similar to that observed macro- 
scopically purely from first-principles transport simula- 
tion and a semi-classical model with simple assumptions. 
Moreover, Rossnagel and Kuan find that resistivity de- 
creases back to roug hly the Cu-only value upon oxidation 
of the interface [lj] . Thus the elimination of the Ta con- 
ducting path actually improves conduction, just as we 
would expect. 

However, as noted, our simple analysis gives only 



FIG. 2: Increasing Cu sheet resistance as a function of Ta 
overlayer thickness, as measured in Ref. [13 ]. Also shown are 
fits to Eq. 0, with Tcu at our calculated value of 0.22 as 
well as a better fit value of 0.54. See discussion in the text 
regarding this quantitative discrepancy. 



4% resistivity augmentation for this geometry, com- 



pared with 10% observed by Rossnagel and Kuan [14 1. 
There are three possible reasons for the quantitative dis- 
agreement on the magnitude of resistivity augmentation. 
First, there may be an underestimate in this work of the 
transmission probability across the interface. Resistivity 
augmentation of 10% for a single Ta overlayer would re- 
quire an effective 1 — p ss 0.3, meaning that transmission 
probability should be about 54%. Note, though, that 
such high T is borderline inconsistent with our calcula- 
tion Mcu/Mra — 1-8. This could be explained by an 
increase in the Fermi level in Ta, or by the small crys- 
tallite size in Ta, leading to many interfaces with other 
orientations of Ta for which M is different. 

Using pi ~ p s may well underestimate the probability 
of diffuse reflection at the interface. Then our calculation 
is a lower limit on scattering due to the liner layers. The 
Cu/Ta epitaxy in real samples introduces non-specular 
reflection from the Ta to within the angles allowed by 
the epitaxial periodicity, even for a perfectly periodic in- 
terface. Finally, it is possible that the electrons would be 
subject to non-specular reflection due to disorder in the 
interface. These last two explanations are also consistent 
with the thickness dependence in Fig. [2 because for Ta 
thickness less than 1 nm, the interface is not fully estab- 
lished. A more detailed study of angular information on 
the reflected electrons at the interface would tell us the 
severity of non-specular reflection at the interface, and 
therefore clarify the last two possibilities. 

In summary, we have calculated transmission proba- 
bility across Cu/Ta liner layer interfaces and developed 
a simple model that matches qualitatively the behavior 
of resistivity. This work suggests that liner layers signif- 
icantly degrade the conductivity of nanoscale wires and 
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that the use of insulating liner layers to reduce the trans- 
mission probability could significantly improve intercon- 
nect conductivity 0]. 
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